Orbital order and magnetism of FeNCN 
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Based on density functional calculations, we report on the orbital order and microscopic magnetic 
model of FeNCN, a prototype compound for orbital-only models. Despite having a similar energy 
scale, the spin and orbital degrees of freedom in FeNCN are only weakly coupled. The ground-state 
configuration features the doubly occupied d 3z 2_ r 2 (ai g ) orbital and four singly-occupied d orbitals 
resulting in the spin S = 2 on the Fe +2 atoms, whereas alternative (E' g ) configurations are about 
75 meV/f.u. higher in energy. Calculated exchange couplings and band gap are in good agreement 
with the available experimental data. Experimental effects arising from possible orbital excitations 
are discussed. 
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The relationship between spin and orbital degrees of 
freedom in transition-metal compounds is well estab- 
lished on both phenomcnological and microscopic levels. 
Phenomenologically, the orbital states are described by 
pseudospin operators and can be treated using diverse 
techniques developed for spin Hamiltonians [lj. Micro- 
scopically, the orbital pattern determines the superex- 
change couplings that are the main driving force of mag- 
netism in insulators. The opposite effect, the influence of 
magnetism on the orbital state, is less universal [2] than 
initially expected Particularly, recent computational 
studies of model orbitally-ordered materials Q indicated 
the important role of lattice distortions in stabilizing spe- 
cific orbital states. Although spin and orbital degrees of 
freedom are inherently tangled, there has been consid- 
erable theoretical effort in exploring orbital-only models 
with no spin variables involved 4 6]. In this paper, we 
will present a compound that features intrinsically weak 
coupling between spins and orbitals, and may be a feasi- 
ble experimental probe for such orbital-only models. 

As a model compound, we consider the recently discov- 
ered iron carbodiimide FeNCN containing Fe +2 cations 
that form layers of close-packed FeN6 octahedra in the 
ab plane Q. The layers are connected via linear NCN 
units (Fig. [T]). Experimental information on FeNCN is 
rather scarce. The compound is a colored (dark-red [7| or 
brown Q ) antiferromagnetic insulator with the Neel tem- 
perature of 345 K [7| . The electronic structure of FeNCN 
was studied by Xiang et al. |8[ who arrived at a puzzling 
conclusion on the dramatic failure of conventional density 
functional theory (DFT)+C7 methods that were unable to 
reproduce the experimentally observed insulating ground 
state. We will show that this failure is caused by a sub- 
tle effect of competing orbital states. Such orbital states 
are readily elucidated in a careful DFT-based study, and 
reveal an unusually weak coupling to the magnetism. 

Our DFT calculations are performed in a full-potential 
code with a local-orbital basis set (FPLO) [9]. We used 
the experimental crystal structure, the local density ap- 
proximation (LDA) with the exchange-correlation po- 




FIG. 1. (Color online) Left panel: crystal structure of FeNCN. 
Right panel: close-packed layer of FeNg octahedra (top), the 
single octahedron squeezed along the three-fold c axis (bottom 
left), and a sketch of the ground-state d 6 electronic configu- 
ration of Fe +2 (bottom right). 



tential by Perdew and Wang [lCf, and well-converged k 
meshes of 1500 — 2000 points in the symmetry- irreducible 
parts of the first Brillouin zone. The application of 
a generalized-gradient-approximation (CCA) exchange- 
correlation potential led to quantitatively similar results. 

The LDA energy spectrum for FeNCN (Fig. [2]) strongly 
resembles that of iron oxides. Nitrogen 2p states form va- 
lence bands below —2 eV, whereas Fe 3d states are found 
in the vicinity of the Fermi level. The energy spectrum is 
metallic due to the severe underestimation of electronic 
correlations in LDA. 

The local picture of the electronic structure stems from 
the crystal-field levels of Fe +2 with the electronic config- 
uration d 6 . The octahedral local environment induces 
the conventional splitting of five d states into ti g and 
e g levels. This primary effect is accompanied by a weak 
trigonal distortion that further splits the t2 g states into 
the aig singlet and e' g doublet (Fig. [IJ bottom). The 
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FIG. 2. (Color online) Top: LDA density of states for FeNCN. 
Bottom: Band structure (thin light lines) and the fit with the 
tight-binding model (thick dark lines). The Fermi level is at 
zero energy. 



balance between the a\ g and e' g states is determined by 
fine features of the local environment. According to sim- 
ple electrostatic arguments, the squeezing of the octahe- 
dron along the three-fold axis should slightly favor the a\ g 
state (note the N-Fe-N angles a = 95.9° and /3 = 84.1° 
in Fig. Q]), which is represented by d3 Z 2_ r 2 orbital in the 
conventional coordinate system (z along the c axis). 

To quantify the orbital energies, we fit the Fe 3d 
bands with a tight-binding model based on Wannier func- 
tions adapted to specific orbital symmetries (Fig. [3]). 
The fit yields e aig = -0.29 eV, e dg = -0.30 eV, and 
e e = 0.58 eV. The ti g — e g splitting of about 0.9 eV is 
typical for 3d systems, whereas the energy separation of 
10 meV between the a\ g and e' g states is very small and 
opposite to the naive crystal-field picture. The difference 
may arise from covalency effects and/or long-range in- 
teractions inherent to solids (note a similar example in 



Ref. [llj) 



Similar to iron oxides, FeNCN is expected to feature 
the high-spin state of Fe +2 (the non-magnetic low-spin 
state apparently contradicts the experimental magnetic 
response reported in Q). Therefore, in the local pic- 
ture five out of six d electrons occupy each of the d or- 
bitals, whereas the sixth electron takes any of the a\ g or 
e' g orbitals, thereby creating orbital degrees of freedom. 
Since the Mott-insulating state implies integer orbital oc- 
cupations, the ground state of FeNCN should feature one 
doubly-occupied and four singly-occupied orbitals. The 
nature of the doubly occupied orbital is determined by 
the energy difference between a\ g and e' g and, more im- 
portantly, by correlation effects. 
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FIG. 3. (Color online) Wannier functions based on the a 



and e' g orbitals. 
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To account for correlation effects in FeNCN, we use 
the DFT+Z7 method that treats strong electronic corre- 
lations in a mean-field approximation valid for insulators. 
In contrast to Ref. [8| reporting the half-metallic solution, 
we readily obtained insulating solutions by explicitly con- 
sidering orbitally-ordered configurations. The failure of 
the previous computational work is likely related to spu- 
rious solutions arising from random starting configura- 
tions. Such solutions are not converged in charge due 
to the charge shuffling effect in a (half)-metal. To over- 
come this problem, we first performed calculations with 
a fixed occupation matrix (i.e., fixed the orbital configu- 
ration) and later released this matrix to allow for a fully 
self-consistent procedure. A similar approach has been 
used in previous computational studies [ill [l2[ , and was 
shown to be vital for the proper treatment of systems 
with orbital degrees of freedom. 

The input parameters of the DFT+U method, the on- 
site Coulomb repulsion (Ud) and exchange (Jd), are eval- 
uated in a constrained LDA procedure [13j implemented 
in the TB-LMTO-ASA code Q. We find [/Lmto = 
6.9 eV and Jd = 0.9 eV, whereas a comparative calcula- 
tion for FeO yields similar values of U^ MT =7.1 eV and 
Jd = 0.9 eV. The magnitude of electronic correlations in 
FeNCN is, therefore, the same as in Fe +2 oxides. By 
contrast, a somewhat reduced Ud parameter was found 
in CuNCN (6.6 eV vs. 9 - 10 eV in Cu+ 2 oxides) and as- 
cribed to sizable hybridization between the Cu 3d and N 
2p states [Hj]. In FeNCN, such a hybridization is rather 
weak (Fig. [2]). 

Since the two Fe sites in the unit cell of FeNCN are 
crystalographically equivalent, we restrict ourselves to 
ferro-type orbital configurations featuring the same dou- 
bly occupied orbital on both Fe sites. Starting from dif- 
ferent occupation matrices, we were able to stabilize sev- 
eral solutions [3]. The ground-state configuration fea- 
tures the doubly occupied a\ g orbital and is further re- 
ferred as Ai g . The E' g configurations with two electrons 
on either of the e' g orbitals lie higher in energy for about 
75 meV/f.u. This energy difference is nearly independent 
of the specific Ud value. It is also possible to put two elec- 
trons on one of the e g orbitals, but such configurations 
are highly unfavorable (0.81 eV/f.u. above the ground 
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TABLE I. Exchange couplings (in K) calculated for the 
ground-state (Ai g ) and one of the higher- lying (E' g ) configu- 
rations. Interatomic distances are given in A (see also Fig. [TJ . 
The on-site Coulomb repulsion parameter is Ud — 7 eV. 
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state) in agreement with the large <2g — e g splitting of 
0.9 eV in LDA. The artificial low-spin configuration with 
all six electrons on the t2 g orbitals has an even higher 
energy of 3.6 eV/f.u. above the A\ g ground state. 

The lowest-energy spin configurations are antiferro- 
magnetic, irrespective of the A\ g or E' g orbital state. The 
energy spectra are quite similar [l7j |. although the band 
gap (-Bgap) for the A\ g orbital configuration is systemat- 



ically higher than for any of the E' configurations: for 
example, at U d = 7 eV £ gap = 2.95 eV and 2.30-2.50 eV 
for A\ g and E', respectively. The change in Ud causes a 
systematic shift of the band gaps [18|. While the lack 
of the experimental optical data prevents us from tuning 
the Ud parameter against the experimental band gap, we 
note that the calculated E gap values agree well with the 
dark-red color of FeNCN. 

We now investigate the interplay between spin and or- 
bital degrees of freedom in FeNCN. To evaluate magnetic 
couplings, we doubled the unit cell in the ab plane, and 
calculated total energies for several spin configurations. 
These total energies were further mapped onto the clas- 
sical Heisenberg model yielding individual exchange in- 
tegrals Ji. We evaluated the nearest-neighbor coupling 
J a i, in the ab plane as well as the nearest-neighbor and 
next-nearest-neighbor interplane couplings J c and J' c , re- 
spectively (Fig. [TJ . Further couplings are expected to be 
weak due to negligible long-range hoppings in our tight- 
binding model. 

Surprisingly, the calculated exchange couplings listed 
in Table [TJ depend only weakly on the orbital order. 
Both Ai g and E' g orbital configurations induce lead- 
ing antifcrromagnetic (AFM) exchange J c via the NCN 
groups. The intraplane coupling J ao is ferromagnetic 
(FM), whereas J' c is AFM. The resulting spin lattice is 
non-frustrated, and features AFM long-range order with 
parallel spins in the ab plane and antiparallel spins in 
the neighboring planes. This prediction awaits its verifi- 
cation by a neutron scattering experiment. To test our 
microscopic magnetic model against the available exper- 
imental data, we simulated the magnetic susceptibility 
using the quantum Monte-Carlo loop algorithm [19| im- 
plemented in the ALPS package 120]. The susceptibility 
was calculated for a three-dimensional LxLxL finite lat- 
tice with periodic boundary conditions and L = 12 21 1. 
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FIG. 4. (Color online) Experimental magnetic susceptibility 
of FeNCN 0] and the simulated curve for J c = 46.5 K, J ab = 
— 16.3 K, and J' c = 2.3 K. The deviations at low temperatures 
are due to an impurity contribution and/or anisotropy effects. 



The simulated magnetic susceptibility was compared 
to the experimental data from Ref. (Fig. 13). The data 
above 100 K and, particularly, the transition anomaly at 
Tm = 345 K are perfectly reproduced with J c — 46.5 K, 
J a b = —16.3 K, and J' c = 2.3 K in remarkable agreement 
with the calculated exchange couplings for the Ai g or- 
bital configuration (Table UJ) [22]. The low-temperature 
upturn of the experimental susceptibility violates the be- 
havior expected for a Heisenberg antiferromagnet, and 
signifies an impurity contribution or effects of exchange 
anisotropy that are not considered in our minimum mi- 
croscopic model. 

The three-dimensional magnetism of FeNCN is some- 
what unexpected considering the seemingly layered na- 
ture of the crystal structure (Fig. [TJ). The leading ex- 
change is antiferromagnetic and runs between the lay- 
ers, although the nearest-neighbor interlayer distance of 
4.70 A is much longer than the intralayer distance of 
3.27 A. The unusually strong interlayer exchange orig- 
inates from the peculiar nature of the NCN units that 
feature a strong 7r-bonding and mediate hoppings be- 
tween the neighboring layers. This effect has been illus- 
trated by sizable contributions of both nearest-neighbor 
and second-neighbor nitrogen atoms to Wannier func- 
tions in CuNCN [15| . Similar contributions are found for 
the e g orbitals in FeNCN. 

The intralayer interaction is a conventional combina- 
tion of the direct exchange and Fe-N-Fe superexchange 
that result in a weakly ferromagnetic coupling (the Fe- 
N-Fe angles are 95.9°, i.e., close to 90°). The spin lattice 
of FeNCN reminds of another transition-metal carbodi- 
imide, CuNCN, where the long-range antiferromagnetic 
superexchange mediated by the NCN groups was also re- 
ported [Hj|. The difference between the two compounds 
is the strong Jahn- Teller distortion in CuNCN that splits 
the close-packed layers of transition-metal octahedra into 
structural chains running along a, with the ferromag- 
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netic coupling resembling J a b in FeNCN. By contrast, 
FeNCN is not subjected to a Jahn- Teller distortion, and 
the ground-state orbital configuration is largely stabilized 
by electronic correlations. Indeed, there is no clear pref- 
erence for a certain orbital state on the LDA level. 

According to our results, the energy scales for the spin 
and orbital degrees of freedom in FeNCN are compara- 
ble, about 52 meV/f.u. and 75 meV/f.u., respectively. 
However, magnetic couplings weakly depend on the or- 
bital configuration keeping spins and orbitals nearly de- 
coupled. This effect is explained by different d states 
responsible for the orbital and magnetic effects. The or- 
bital degrees of freedom are operative in the ti g subspace, 
whereas magnetic couplings are largely determined by 
the e g orbitals featuring stronger intersite hoppings. The 
decoupling of spin and orbital variables along with the 
low energy scale of the competing orbital states suggest 
that the orbital-only physics can be probed in FeNCN. 

The characteristic energy scale of 75 meV/f.u. corre- 
sponds to temperatures around 850 K, and implies that 
the E' g orbital states may emerge at elevated tempera- 
tures. Although FeNCN, alike all transition-metal car- 
bodiimides, is thermodynamically unstable 
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be maintained up to at least 680 K, which is the prepa- 
ration temperature reported in Ref . 0] ■ Other options of 
activating the E' g orbital states of FeNCN are the applica- 
tion of pressure and laser irradiation. The latter has been 
successfully used for melting the orbital order in several 
prototype orbital systems [24J and could be applied to 
FeNCN as well. If the switching of the orbital state is 
possible, the anticipated effect is a sizable reduction in 
the band gap (for at least 0.4 eV), while the structure will 
probably adjust to the new orbital state. However, it is 
more likely that several competing E' g states will form an 
orbital liquid, thereby maintaining the high symmetry of 
the crystal structure. To probe such effects, further ex- 
perimental work on FeNCN is highly desirable. We also 
mention an isostructural compound CoNCN |25|, where 
orbital degrees of freedom arising from Co +2 (d 7 ) are ex- 
pected. 

In summary, we have shown that FeNCN presents an 
unusual example of weakly coupled spin and orbital de- 
grees of freedom acting on a similar energy scale. The 
ground-state orbital configuration features two electrons 
on the a\g orbital, in agreement with the simple electro- 
static arguments, but in contrast to the LDA-based ex- 
pectations. The calculated properties, such as the band 
gap, exchange couplings, and Neel temperature, are in 
very good agreement with the experiment. We have also 
remedied the failure of the recent computational work [8j 
and confirmed the remarkable performance of DFT+J7 
techniques applied to Mott insulators with orbital de- 
grees of freedom. 

We are grateful to Deepa Kasinathan and Oleg Jan- 
son for fruitful discussions. We also acknowledge Richard 
Dronskowski and Andrey Tchougreeff for drawing our at- 
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FIG. 5. Atomic- and orbital-resolved density of states for the A\ g (left panel) and one of the E' g (right panel) orbital 
configurations of FeNCN. The Fermi level is at zero energy. The on-site Coulomb repulsion parameter is U& = 7 eV. 



